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1 Introduction 

Despite recent trends towards ever larger computers, I believe that accurate 
and reliable lattice simulations of QCD are possible on quite modest computers, 
possibly even on the workstation or PC sitting on your desk. If true, this 
will have a profound effect on the way we deal with strong interactions, both 
experimentally and theoretically. In these lectures I outline the reasons behind 
my belief. 

The central issue is the reliability of lattice simulations of QCD. Lattice 
simulations, like experiments, have two sources of error: statistical errors due 
to our use of Monte Carlo methods for doing the simulations; and systematic 
errors due to the various approximations we make to simplify the theory to 
the point that it is numerically tractable. Statistical errors are pretty much 
understood, and under control; I do not discuss them further in these lectures. 
Instead, I concentrate on systematic errors. It is the control of systematic errors 
that determines the scale and reliability of contemporary QCD simulations. 

I begin in Section 2 with a discussion of the main sources of systematic error 
in lattice simulations. Of these errors those due to the finite lattice spacing are 
the most important. I describe a general strategy for systematically improv- 
ing lattice lagrangians to remove lattice-spacing errors. This is an alternative 
to simply reducing the lattice spacing that promises to be much more efficient. 
This strategy relies upon weak-coupling perturbation theory, and so in Section 3 
I discuss the extent to which perturbation theory can be trusted in lattice QCD. 
Contrary to much conventional lore lattice perturbation theory seems to work 
well even at distances as large as 1/2 fm. In Section 4 I discuss another ob- 
stacle to the perturbative improvement of lattice QCD. The problem is with 
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"tadpole" terms that arise in perturbation theory when computing corrections 
to the QCD lagrangian. The tadpoles are usually too large to be dealt with 
perturbatively. All such contributions have a common origin, and I describe a 
simple nonperturbative procedure for dealing with them. In Section 5 I illus- 
trate all of these techniques with recent simulation results for the T family of 
mesons — an area where lattice QCD has already largely succeeded. Finally, in 
Section 6, I summarize the central result of these lectures. 

2 Systematic Errors in Lattice QCD 

There are three important sources of systematic error in lattice QCD: 

a) finite volume errors that arise because simulations cover only a finite vol- 
ume V = L A in space-time; 

b) errors due to extrapolations from large quark masses, where the simu- 
lations are easy, down to realistic quark masses for u and d quarks, or, 
cquivalcntly, from large m 2 oc ?>?„.,; to m 2 ^ = (140 MeV) 2 ; 

c) errors due to the finite lattice spacing a in the grid used to approximate 
continuous space-time. 

The central question in numerical QCD then is: How large must L be, and how 
small m n and a, for simulations to be realistic (and useful)? The answer to this 
question determines the cost and feasibility of lattice simulations. 

The cost of a lattice simulation is at least proportional to the number of 
sites in the grid, 



Usually, however, it is more costly because of a phenomenon called "critical 
slowing down" that tends to slow the convergence of the iterative algorithms 
used in lattice simulations. A more realistic estimate of the cost is 



where typically a and (3 are somewhere in the range 0-1 (usually closer to 1 
than 0). Obviously this cost is highly sensitive to the values of a, m w , and L; 
for example, in the worst case, the cost is proportional to the seventh power of 
1/a. Thus it is critically important to make a and m T as large as possible, and 
L as small as possible. We cannot afford to be sloppy or to "play it safe." 

Here I discuss each of these systematic errors in turn. However my emphasis 
is on the finite-a errors as these are by far the most crucial determinants of the 
cost (see Eq. (0)), and they are also the most misunderstood of the systematic 
errors. 
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2.1 Finite L 



Clearly when we simulate a hadron, like the pion or the proton, that is roughly 
2 fm in diameter, we must have L at least as large as 2 fm.The nature of 
the systematic error when L > 2 fm is easily understood. Lattice simulations 
generally use periodic boundary conditions to define the differential operators in 
the action. This means that one is never studying just a single hadron on a grid, 
but rather an infinite crystal of hadrons consisting of the hadron and images 
of it, induced by the boundary conditions, at intervals of L in each direction. 
The properties of our hadron are affected by interactions with all of its images. 
Its energy, for example, is shifted by the interaction energy between the hadron 
and its images. The interaction energy between two hadrons a distance r apart 
is 

K- • (3) 

r 

Thus the finite volume error falls off exponentially quickly with increasing L 
provided L is large enough. This will not be the case when L is so small that 
the hadron overlaps with its image; then one finds powerlaw effects related to 
the geometrical volume of the overlap. But for L > 2 fm we expect 

finite L errors oc e'"^ L « e' L/1 fm (4) 

This formula suggests that finite volume errors for hadrons ought to be 10% 
or less for L's of order 2 fm or larger, and this seems to be the case in current 
simulations. Even if 2 fm is too small, the exponentially-falling errors vanish so 
quickly with increasing L that it is almost certain that 2.5 or 3 fm will suffice. 
The combination of exponentially small errors and powerlaw cost means that 
finite-volume effects are both tractable and fairly well understood. 

2.2 Large m w 

The numerical intractability of small quark masses means that simulations rarely 
have pion masses smaller than about 500 MeV. (Remember that oc m u ^ 
is a measure of the quark mass.) As a consequence, it is standard practice in 
lattice simulations to present results as functions of m„, extrapolating down to 
the physical mass (140 MeV). The reliability of this procedure depends upon 
how sensitive the simulation results are to the pion mass. 

We have, in chiral perturbation theory, a reasonably reliable tool for explor- 
ing the m v dependence of different quantities. With chiral perturbation theory 
we can compute the m T dependence of various quantities. Some, like the nu- 
cleoli mass, are quite sensitive to when « 500 MeV. Others, like the pion 
decay constant /„•, are relatively insensitive. This suggests that extrapolations 
in ra^ ought to work well for / w , while extrapolations of the nucleon mass could 
lead to problems. 
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Chiral perturbation theory can help us decide which quantities are least 
sensitive to m, , and we might wish to concentrate our attention on these. More 
importantly, we should push on to smaller m,'s. Large errors are typically 
be of order 

large errors 2 . (5) 

This suggests that 's of order 250 MeV or less ought to result in systematic 
errors that are less than 10%. Given that the cost grows only as 1/m 2 ., it 
should be possible to work at these lower masses, especially if we are able to 
use significantly larger lattice spacings. 

2.3 Finite a — Classical Field Theory 

In thinking about errors due to the finite lattice spacing, I find it useful to look 

first at a classical field theory on a lattice. Consider, for example, Poisson's 
equation: 

d 2 4>{v)=p(v). (6) 
On a grid, we replace derivatives by differences, for example 

flgtfr) -> ^ + ^)-20(r) + ^(r- a x) ; (?) 

and specify <fi only at the nodes of a lattice with spacing a. 

Obviously a must be small compared with the size of any important structure 
in p or <f>. If this is the case, the finite-a errors fall off like a power of a. This is 
easily demonstrated. The difference operator that replaces the second derivative 
is defined by 

a 2 A^c/)(r) = (j>(r + ax) - 20(r) + 0(r - ox) 

= (c ac% -2 + e~ ad *) 0(r). (8) 

Expanding the exponentials, we find that 

Ai 2 )^(r)«^ + ^ + ...^(r). (9) 
Thus by replacing d 2 with A^ 2 ) we make errors of order 

2 2 

finite a errors « ^ (10) 

where p is the typical wavenumber (ie, momentum) in the Fourier transform 
of (j){x). 

There are two lessons to learn from this example. The first is that only a 
small number of lattice points is needed to obtain precision of order 10%. For 
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example, if <j> consists of a single smooth bump of width w, then a lattice spacing 
a rts w/3 gives roughly 10% accuracy: p « 27r/A, where Ak2id« 6a, implies 
that a 2 p 2 /12 « 0.08. A general rule of thumb is that it takes roughly 3 d grid 
points per d-dimensional bump to begin to model that bump accurately on a 
grid. 

The second lesson is that one should not reduce the lattice spacing much 
below w/3 when trying to improve the accuracy of the grid approximation. This 
is because reducing the lattice spacing greatly increases the cost of the simulation 
while only modestly improving the errors. For example, to reduce 10% errors to 
1% errors using our discretization of the three-dimensional Poisson's equation 
one would have to decrease a by a factor of 3, thereby increasing the number of 
grid points by a factor of 3 3 = 27. Taking account of critical slowing down, the 
cost of the simulation would increase by a factor of 80-200. 

A far more efficient way to reduce errors is to improve the discretization, 
keeping the lattice spacing large. For example, from our discussion above, 

fi 2 A ( 2 ) _ fl2 ^ _ fl4 dx _ 

' ~~ x 12 360 

12 90 v ' 

and this implies that the lattice equation 

E ( A f ) - ^t 1 ) ^ = >W ( 12 ) 

is accurate up to corrections of order (op) 4 /90. This improved lattice equation 
gives roughly 1% accuracy with a — w/3. The added complexity of the equation 
increases the cost of simulating, but only by a factor of 1.5 or 2. 

Thus for classical fields we conclude that coarse lattices are optimal — perhaps 
3 d grid points per bump in the answer. High precision is obtained by improv- 
ing the lattice operators from which the field equations are constructed, not by 
reducing the lattice spacing. 

Exercise: Compare exact derivatives with lattice derivatives for (f>(x) — exp(—x 2 /2) 
in one dimension. Try lattice spacings of a = 0.5 and 1, and study the second 
derivative at x — 0. Compare the errors obtained using the leading-order lattice 
derivative and the improved derivative. Check how these errors scale with a. 
The lattice derivatives can be corrected to all orders in a by continuing the 
process outlined above. Compute the error in the infinite-order lattice approx- 
imation of d 2 cj>(0). In general, for what sorts of function 4>{x) will the error 
be exponentially small, power-law suppressed, or identically zero? (Hint: Use 
Fourier transforms. Fourier transforming the lattice operators is simple. In par- 
ticular, the infinite-order lattice approximation to d 2 has transform — k 2 , where 
k is the wavenumber.) 
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This last exercise shows that our lattice derivatives can never be perfect; typi- 
cally there are residual non-powerlaw differences between them and continuum 
derivatives. Examine the possibility of removing even these errors when r ^> a in 
our lattice Poisson's equation with a source p(r) that is localized around r = 0. 
In particular, consider adding local counterterms to the source: p(r) — > p + 
Co <5(r) + C2 d 2 S(r) + • • •. (Note the analogy with multipole expansions.) 

Exercise: Compute the ground-state energy of a nonrelativistic particle trapped in 
a box of length L using the Shrodinger equation with discretized derivatives of 
the sort discussed above. (This is easily done analytically using Fourier series 
for the eigenfunctions.) Verify the "3 points per bump" rule and explore the 
utility of improving the discretization as we did above. 

As a variation, consider adding a harmonic potential to the hamiltonian and 
solving the eigenvalue problem numerically. Defined on a grid, the hamiltonian 
operator becomes a matrix, and the energy eigenvalue problem a simple matter 
of determining the eigenvalues of that matrix. Again explore the relative merits 
of reducing the lattice spacing versus improving the discretization. Try the 
problem in 2 or more dimensions as well. 

2.4 Finite a — Quantum Field Theory 

The analysis of finite-a errors for quantum theories is complicated by the fact 
that quantum fields are rough at all length scales. This roughness is due to 
quantum fluctuations. In perturbation theory, it is responsible for the ultraviolet 
sensitivity and infinities of loop diagrams. More generally, it seems to call into 
question the validity of a discrete approximation. Indeed it is hard to imagine 
defining continuum derivatives for an infinitely rough field, let alone discrete 
approximations to them. 

Somewhat surprisingly, finite-a errors in quantum theories arc rather similar 
to those in classical theories. This is because only the long- wavelength structure 
is physical on a lattice. In general, any long-wavelength, low-momentum probe 
is sensitive only to fields averaged over a region of order the size of the probe. 
This averaging smoothes out quantum fluctuations. Consequently the long- 
wavelength or infrared behavior of a quantum theory is insensitive to the details 
of its short- wavelength or ultraviolet behavior; and thus there are infinitely many 
theories, each with different ultraviolet dynamics, that give identical or nearly 
identical infrared physics. We take advantage of this fact when we construct 
our lattice theory: we replace the correct ultraviolet dynamics of the continuum 
theory with the incorrect, but numerically tractable, dynamics of the lattice. 
By changing its bare coupling constants, we are able to adjust the dynamics of 
the lattice theory so that its infrared behavior is the same as in continuum up 
to corrections that vanish with the lattice spacing a — just as in classical lattice 
theories. 
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To see how this works consider a <fr field theory defined on a lattice: 

£ lat = -i]>>A( 2 )0+imV + -lA0 4 (13) 

This lagrangian, used in a path integral, defines a euclidean lattice approxima- 
tion to the continuum quantum theory. Classically, large distortions due to the 
grid arc restricted to the ultraviolet modes; but, in the quantum theory, the in- 
frared sector is strongly affected as well since it couples to the ultraviolet modes 
through quantum fluctuations (ie, loops in perturbation theory). For example, 
distortion of the ultraviolet modes changes the mass renormalization, thereby 
shifting the renormalized mass of the <j> particle — a significant modification of 
the infrared behavior of the theory. 

Introducing a grid causes 0(1) changes in the behavior of the quantum 
theory at all length scales, but the continuum behavior at large distances can 
always be restored by shifting or "renormalizing" the bare coupling constants. 
The renormalized mass of the <j> particle, for example, is easily restored to its 
proper value by adjusting the bare mass m in the lattice action. The grid is 
nothing more than an ultraviolet cutoff, restricting the path integral to momenta 
smaller than ir/a. Renormalization theory tells us that the errors caused in, say, 
a scattering amplitude by introducing a finite cutoff A are suppressed by powers 
of Pext/A, where p cv ± is the largest external momentum in the amplitude. This 
is true provided the coupling constants of the theory are shifted appropriately 
from their continuum values. Thus, for our lattice </> 4 theory, there are particular 
values of m and A such that the lattice result for any physical quantity is the 
same as in the continuum theory up to corrections that vanish as powers of ap, 
where p is a momentum associated with the classical (as opposed to quantum) 
size and scales relevant to that quantity. Generally the values needed for A 
and m depend upon the lattice spacing a (the couplings are said to "run" with 
varying a). 

The problem of dealing with the effects of complex short-distance structure 
on long-distance behavior arises in many contexts, both quantum and classical. 
A classical example is the formulation of electrostatics inside a dielectric medium 
(a piece of glass, for example). The microscopic electric field inside a dielectric 
is very complicated and fluctuates rapidly even over distances as short as 1 A. 
Usually, however, it is not the microscopic field that is physically relevant. Any 
macroscopic probe of the electric field sees only the average field, averaged 
over a large region that includes many atoms. This averaging smoothes out 
the rapid fluctuations, producing a relatively smooth macroscopic field that 
is described by simple equations. The macroscopic equations are universal in 
form; the type and nature of the atoms making up the dielectric are largely 
irrelevant. The only information needed about the microscopic structure of 
the material is contained in the numerical value of its dielectric constant. The 
dielectric constant is analogous to the bare couplings in our lattice theory; and 
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the standard formalism for dielectrics comes from a straightforward application 
of renormalization ideas to the classical problem. 

Despite quantum fluctuations, the finitc-a errors in our quantum theory are 
quite similar to those in the classical theory. Also as in the classical case, we can 
correct the theory to systematically reduce these errors. Renormalization theory 
tells us that all errors of order {ap) n can be removed by adding local interac- 
tions of dimension 4 + n and lower to the lattice lagrangian with appropriate 
couplings (and readjusting any couplings already present in the lagrangian). 
The correction terms are local because they are correcting for the physics at 
distances shorter than the lattice spacing. Only terms that preserve the sym- 
metries of the theory need be added, and so there aren't very many for small n. 
Our lattice <fi A theory, for example, has no possible correction term of dimension 
five. This means that the leading finitc-a errors are quadratic in a. There arc 
only two dimension six correction terms, and so the 0(a 2 p 2 ) errors in £\ at are 
removed by adding 

^^£^ 2) )^+^ 6 (14) 
and by further shifting m(a) and A(a). 

Exercise: There are actually two other dimension-6 operators that are consistent 
with the symmetries of the lattice (f> 4 theory: (A' 2 ^) and 3 A^ 2) 0- We do 
not include these in our corrected lagrangian because they are redundant. These 
operators can always be removed from the lagrangian by making the replacement 

4>{x) -» <t>(x) = 4>{x) + a 2 h 3 (z) + ah A (2) <Kz) (15) 

in the lagrangian, where /i and /2 are suitably chosen functions of the bare mass 
and couplings. The path integral for the lattice theory is an ordinary multidi- 
mensional integral, and this field redefinition is merely a change of integration 
variable; it cannot alter the physical content of the theory. 

Show that /i and fi can always be adjusted so that both redundant operators 
cancel when <f> is transformed in the lagrangian. This is easy to see at tree- level 
(ie, in the classical theory). What complications arise in one- loop order? For 
example, what happens to the jacobian associated with the field transformation? 
While a field transformation of this sort has no effect on physical quantities, 
it usually changes off-shell Green's functions. Thus the Green's functions of 
our improved theory may still deviate from the continuum Green's functions in 
order a 2 . However particle masses, on-shell scattering amplitudes, and the like 
will all be accurate through order a 4 . This is all that we need. Examine the 
differences between on-shell and off-shell quantities using tree-level perturbation 
theory for equivalent theories with and without the redundant operators. 
An important step in discretizing a quantum field theory for nonperturbative 
analyses is to verify that the lattice version of the theory is stable. This is 
an issue when improving the lagrangian. For example, if it happened that the 
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(f) 6 coupling, d(a), was negative, we would have to worry about stabilizing the 
improved theory against (j> — > — oo. Luckily we have tremendous latitude in the 
design of our corrections. For example, we can freely add 4> s , 4> 10 . . . terms to the 
lagrangian without affecting simulation results through order (ap) 4 (provided all 
couplings are readjusted appropriately). Such terms could be used to stabilize 
the discrete theory. The presence of an instability in a lattice theory usually 
becomes obvious when one tries to simulate the theory numerically. 

We conclude that there exists a set of values for the couplings m(a), A(a), 
c(a), and d(a) such that the lattice theory gives the same physical results as 
the continuum theory up to corrections of order (ap) 4 when p <C ir/a. In the 
classical limit of our lattice theory, only the c(a) coupling is necessary: taking 
c(a) = 1 cancels the 0(a?) errors from the derivative operator in the original 
lagrangian. In the quantum theory c(a) must be shifted, and a nonzero d(a) 
added to cancel 0(a 2 ) errors induced by quantum fluctuations. This raises the 
problem of how we should determine the correct values for these couplings. In 
principle we might run simulations for many different values of the couplings, 
and numerically search for the set that gives correct physics. But a defect of our 
lattice theory is that it has many more coupling constants than the continuum 
theory — usually far too many to tune numerically. This is a serious obstacle to 
the use of improved lagrangians. In many theories, including QCD, the situation 
is salvaged by using perturbation theory to express the extra couplings in terms 
of the original two couplings, m(a) and A(a). In perturbation theory, we expect 



The radiative corrections for such couplings involve only momenta of order ir/a 
and larger since these are the momenta where the lattice theory needs correcting. 
In asymptotically free theories, like QCD, the fundamental coupling constant 
is small and perturbation theory is valid when the lattice spacing is small. In 
such situations perturbation theory can be used to correct the lagrangian, even 
if the long-distance behavior of the theory is highly nonperturbative. 

To summarize, finite-a errors in quantum lattice theories are very similar 
to those for classical lattice theories. In both cases the errors can be removed, 
order by order in a, by correcting the lattice operators. The main difference in 
the quantum case is that the coefficients of the correction terms in expressions 
like Eq. (|ll| ) are renormalized by a-dependent quantum loop effects that are 
very specific to the particular theory and context in which the operator is used. 

Exercise: When the lattice couplings are perturbative, they may be computed by 
matching perturbative results from the lattice theory against the corresponding 
results generated with the continuum theory. If one has n couplings to compute, 
one chooses n low-momentum physical quantities to match. Each quantity is 
computed, order by order in perturbation theory, both in the lattice theory 
and in the continuum theory. The lattice couplings are then adjusted so that 



c(a) 
d(a) 



1 + c 2 (am) A(a) 2 + ■ • • , 
d 3 (am) A(a) 3 H . 



(16) 



the lattice results agree with the continuum results. There are no ultraviolet 
infinities in the lattice calculations, since the theory has a finite cutoff, and 
therefore the bare couplings are all finite numbers. 

The coupling d{a) is easily computed by examining the one-loop scattering am- 
plitudes for 00 — > 0000; and the coupling c(a) can be obtained from the disper- 
sion relation, relating the energy and momentum of a particle, that is obtained 
by locating the poles in the 0^0 propagator. Sketch these calculations to one- 
loop order. 

2.5 QCD on a Lattice 

Our discussion in the previous section indicates that there are really two re- 
strictions on the maximum lattice spacing a that is useful for numerical work 
in QCD: 

1) a must be smaller than any important scale in the hadronic system under 
study. Since light hadrons have diameters between 1.5 and 2 fm, the 
lattice spacing should probably be no larger than about 0.5 fm. Such 
a lattice spacing should suffice for calculations of static properties, like 
masses, magnetic moments, charge radii, etc. Smaller lattice spacings are 
needed if the hadrons have nonzero momenta. 

2) a must be sufficiently small that perturbation theory works well at dis- 
tances smaller than a (or momenta larger than n/a). If this is the case 
we can use perturbation theory to compute the correction terms needed 
in the lagrangian to obtain high precision from coarse lattices; in effect, 
we can fill in the cracks in the lattice. 

The second of these restrictions is the more controversial. There is much con- 
ventional wisdom suggesting that lattice spacings as small as 0.1-0.05 fm might 
be required to satisfy this condition. I will argue, in subsequent sections, that 
o's as large as 0.5 fm are probably alright. This means that lattices can be much 
coarser than is usual today. Before pursuing this crucial point, we must first 
discuss the formulation of QCD on a lattice. 
The continuum lagrangian for QCD is 

C = -^TrF%, (17) 

where 

Ffiu = d„A v - d v A^ + %g[A^A v ] (18) 

is the field tensor, a traceless 3x3 hermitian matrix. The defining characteristic 
of the theory is its invariance with respect to gauge transformations where 

F„„ -f Sl(x) F„ v nix)' 1 (19) 

and Q(x) is an arbitrary x-dependent SU 3 matrix. 
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The standard discretization of this theory seems perverse at first sight. 
Rather than specifying the gauge field by the values of A^x) at the nodes 
of the lattice, the field is specified by variables on the "links" joining the nodes. 
In the classical theory, the "link variable" on the link joining a node at x to one 
at x + afi is determined by the line integral of along the link: 

U^x) = P exp (-ig jf + " A ■ d^J (20) 

where the P-operator path-orders the A^s along the integration path. We 
use t/ju's in place of A^s on the lattice, because it is impossible to formulate 
a lattice version of QCD directly in terms of A^'s that has exact gauge in- 
variance. The U M 's, on the other hand, transform very simply under a gauge 
transformation: 

U^^flix + a^U^flix)- 1 . (21) 

This makes it easy to build a discrete theory with exact gauge invariance. 

Exercise: Show that a gauge invariant object can be associated with any closed path 
on the lattice by forming the product of U^'s and U^'s associated with the links 
of the path, and taking a trace. In general U ti (x) is associated with a link leaving 
site x in (positive) direction /u, while U^(x) is associated with a link entering 
site x from direction fi. The trace of a product of link variables along a closed 
path is known as a "Wilson loop." 

You might wonder why we go to so much trouble to preserve gauge invariance 
when we quite willing give up Lorentz invariance, rotation invariance, etc. The 
reason is quite practical. With gauge invariance, the quark-gluon, three-gluon, 
and four-gluon couplings in QCD are all equal, and the bare gluon mass is zero. 
Without gauge invariance, each of these couplings must be tuned independently 
and a gluon mass introduced if one is to recover QCD. Tuning this many pa- 
rameters in a numerical simulation is very expensive. This is not much of a 
problem in the classical theory, where approximate gauge invariance keeps the 
couplings approximately equal; but it is serious in the quantum theory because 
quantum fluctuations (loop-effects) renormalize the various couplings differently 
in the absence of exact gauge invariance. So while it is quite possible to for- 
mulate lattice QCD directly in terms of A^s, the resulting theory would have 
only approximate gauge invariance, and thus would be prohibitively expensive 
to simulate. Symmetries like Lorentz invariance can be given up with little cost 
because the symmetries of the lattice, though far less restrictive, are still suf- 
ficient to prevent the introduction of new interactions with new couplings (at 
least to lowest order in a). 

We must now build a lattice lagrangian from the link operators. We require 
that the lagrangian be gauge invariant, local, and symmetric with respect to 
axis interchanges (which is all that is left of Lorentz invariance). The most local 
nontrivial gauge invariant object one can build from the link operators is the 



11 



"plaquette operator," which involves the product of link variables around the 
smallest square at site x in the fiv plane: 

TrU^x) ee Tr (U„(x)U v {x + a£)t/t(x + a/t + aO)Ul(x)) . (22) 

To see what this object is, consider evaluating the plaquette centered about a 
point xq for a very smooth weak classical field. In this limit, 

Tr Uf^ w 3 (23) 

since 

Ufj, rs c - l9aA » w 1. (24) 

Given that is slowly varying, its value anywhere on the plaquette should 
be accurately specified by its value and derivatives at xo- Thus the corrections 
to Eq. ((2^) should be a polynomial in a with coefficients formed from gauge- 
invariant combinations of A^(xq) and its derivatives. Thus we expect 

TrU£ l r = 3 - Cl a 4 Tr( 5 F^K)) 2 

- c 2 a 6 Tr (gF^(x )(Dl + Dl)gF^(x )) 

+ 0(a s ) (25) 

where c\ and ci are constants, and is a gauge-covariant derivative. The lead- 
ing correction is order a 4 because F^ u is the lowest-dimension gauge-invariant 
combination of derivatives of A^, and it has dimension 4. It is a simple exercise 
to show that ci = 1/2 and C2 = 1/24 (in the classical limit). 

The expansion in Eq. (^H|) is the classical analogue of an operator product 
expansion. It provides the classical relation between the plaquette operator and 
the local gauge-invariant operators of the continuum theory. This relationship 
is preserved in the quantum theory except that: a) the expansion parameters Ci 
are renormalized by quantum fluctuations; and b) additional terms enter at 
order a 6 and higher. So the lattice equivalent of the continuum lagrangian in 
the quantum theory is 

* Re (i Tr U$* - l) ee -\ Tr F%, (26) 
y a Uq 

up to corrections of order a 2 . Here uq accounts for the quantum renormalization 
of the classical relation; it is usually absorbed into the definition of the coupling: 

9ut=gul. (27) 
Thus the lattice action for QCD is 

S = P Re ( 3 Tr U t^) + constant (28) 
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where 

(3=4~- (29) 

#lat 

Notice that the lattice spacing has disappeared from the action. It is customary 
in lattice simulations to use "lattice units" for which a = 1, thereby completely 
removing a from the simulation. The lattice spacing enters only implicitly, 
through the numerical value of the parameter (3; different (3's correspond to 
different a's since the bare coupling constant is a function of the lattice spacing. 

Exercise: Find correction terms for the classical lattice lagrangian that remove the 
order a 2 errors. Consider, for example, adding a term that involves a x 2a 
rectangular Wilson loops. 

Exercise: In the classical lattice theory there is only one order a 2 correction. When 
quantum effects are included, other operators appear in this order. What opera- 
tors are possible? Remember that these operators must be local, gauge invariant, 
and respect lattice symmetries. 



3 Perturbation Theory: Where are the Limits 
in Lattice QCD? 

As we discussed in the previous section, the cost of numerical simulations of 
QCD is dramatically reduced as the lattice spacing a is increased. The resulting 
loss in accuracy can be minimized by improving the lattice lagrangian, using 
perturbation theory to compute the corrections. Such perturbative improvement 
is only possible if perturbation theory is applicable at momenta of order n/a. 
Thus the critical question as we increase a is whether or not perturbation theory 
works at momenta of n/a and larger. 

Empirical evidence suggests that perturbative QCD applied to continuum 
quantities works down to momenta as low as 1 GeV; for example, t decay pro- 
vides one of the most accurate determinations of the strong coupling constant. 
Thus lattice spacings as large as 0.5 fm ought to work. This result is very 
encouraging, but it is important to confirm it with detailed numerical tests 
using lattice simulations. Indeed, given the tremendous control we have over 
simulations, it is far easier to test perturbative QCD against nonpcrturbative 
simulations than against experimental data. 

One common procedure for testing perturbation theory using simulations is 
to compute the (3 function — a perturbative quantity — by comparing values for 
some physical quantity, like the glueball mass or the string tension, computed at 
various values of the bare coupling. This gives a nonperturbative determination 
of the (3 function to compare with the perturbative prediction. In my opinion, 
this is a bad way to test perturbation theory since failure in such a test could 
reflect one of two possible problems: either perturbation theory is not working, 
or there are large 0(a, a 2 ) errors in the calculation of the physical quantity used 
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to determine the j3 function. It is fairly difficult to rule out the second option 
numerically, and so a negative result is ambiguous. 

A better procedure is to follow phenomenologists, who identify short-distance 
quantities that can be both measured nonperturbatively and calculated in per- 
turbation theory. As lattice theorists, we have a much simpler job than the 
phenomenologist since we can easily design dozens of quantities that are very 
ultraviolet and that are easy to measure in simulations; we don't have to worry 
about experimental cuts, Monte Carlo corrections for hadronization, and all the 
other problems associated with experimental data. The idea is to compute var- 
ious short-distance quantities — for example, vacuum matrix elements of local 
operators like Tr Uj^ q — using numerical simulations and then to compare the 
results with predictions from perturbation theory. 

Careful testing of perturbation theory requires a careful definition of the 
expansion parameter, the strong coupling constant. This is as true in lattice 
calculations as it is in continuum calculations. There are two aspects to the 
problem. First we need to define what we mean by the running coupling con- 
stant a s (q) ("the scheme"), and then we need a procedure for choosing a q 
("the scale") appropriate to the quantity of interest. We deal with these issues 
in the next section. Then, in the following section, we compare perturbative 
results with results from nonperturbative simulations. Note that our choice of 
scheme and scale for a s is completely automatic; there is no fine tuning of the 
perturbative results we use in our comparison with simulations. 

3.1 Defining a s 

Given one definition ai(q) it is easy to create other definitions: for example, 
a 2 (q) = Oii(2q) or a 3 (q) = ai(q) + 3ai(q) 2 + In principle, any of these can 
be used as the expansion parameter for perturbation theory; in practice some 
definitions will be more reliable than others. For example, if ai(q) is a reliable 
choice then cti{q) = ai(q) + 100ai(g) 2 will probably be useless. Expansions in 
terms of either coupling will be the same in lowest order, but higher orders will 
have large negative corrections when 0:4 is used. To minimize such problems, I 
like to use a physical quantity to define the coupling. One such definition that is 
particularly convenient uses the potential between a static quark and antiquark: 



where Cp = 4/3 and q is the momentum transferred between the quark and 
antiquark. I like this definition because there is no ambiguity about the rela- 
tionship between the argument of ay and the momentum scale in the process: 



V(q) EE - 




(30) 



a v (q) = 



the coupling strength of a 
gluon with momentum q. 



(31) 
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This is a useful definition both for continuum and for lattice results. It is easy 
to convert expansions in terms of a^g to ay expansions by usingd 



a- 



MS 



(q) = a v (e 5/6 q){l + 2a v /ir+---}. (32) 



Having settled on ay(q), we now need some way of specifying the scale q. 
The procedure is best understood through an example. Consider the vacuum 
expectation value (1 — | Trl^} computed in Landau gauge. To leading order, 
this quantity is proportional to (Tr^4 2 ) and the leading contribution in lattice 
perturbation theory is 

(l-±Triy LG = ay(q*)2™ 2 / ^ ^ 1 + 0{c? v ) (33) 

where q^ = (2/a) sin(a c^/2), and l/<? 2 comes from a gluon propagator. The 
question is what values for q* makes sense. Given that ay is the coupling for a 
gluon of momentum q, we really want 

,-7r/a ,4 r-x/a >4 

a v {q*) ±1 = ^-a v (q) (34) 

Jn/a T Jrr/a T 

except that the right-hand side diverges at small q. To see why it diverges, we 
can rewrite that side of the equation in terms of ay(/z) for some fixed value 
of fx: 

a ^ {av(n) ~ Poa 2 v (p) Hq 2 / f i 2 ) - PW v {ii) ln 2 (q 2 //i 2 ) -•••}. (35) 

-7T /a Q 

Each term in this expansion separately is finite, but the sum diverges. This is 
not surprising since perturbation theory yields only an asymptotic expansion. 
Furthermore it is inconsistent to keep an infinite number of radiative corrections 
in this expansion while only keeping only one or two orders elsewhere in a 
calculation. So, for present purposes, we should replace the ay's in Eq. (|34|) by 
the first two terms in their expansion in terms of ay (/i). It is then easy to solve 
for q* to obtain: 

/■K la, j4 

ln(<n = " /a / q 2 ■ (36) 

r 7r/a d i q 

IT J (l <f 

Evaluating the integrals one obtains q* = 2.8 /a, and a final prediction of 

(1 -^TtUJlg = 0.97ay(2.8/a) (37) 



1 In practice there is little difference between using o^jg- or using ay, but the latter suggests 
a simple procedure for setting the scale. 
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The analogous quantity in the continuum has a quadratic ultraviolet divergence, 
and so it is not surprising that the average momentum scale q* in this process 
very nearly equals the maximum momentum allowed on the lattice. 

This example illustrates a general procedure. A quantity whose first order 
contribution has the form 

a v (q*) [ ' d 4 qf(q), (38) 

J — tt/ a 

where q is the gluon momentum, has scalei 

a d 4 qf(q) ln(q) 

tt j a 



MO=^ • (39) 

d\f{q) 

n j a 

Each perturbative quantity has its own characteristic scale; the more ultraviolet 
a quantity is, the larger is its q* . And, in an asymptotically free theory, the 
larger q* is, the more accurate is perturbation theory. 



3.2 Testing Perturbation Theory 

To probe the validity of perturbation theory, we need simulation results for 
several short-distance quantities and we need perturbative predictions for those 
quantities. Here we focus on logarithms of the expectation values of planar 
m x n Wilson loops, 

lnW mn = ln(±Tr[/ mxn ), (40) 

because these are easy to simulate and have been analyzed in perturbation 
theory through second order. (We use the logarithms of the loops rather than 
the loops themselves because the logarithms have more convergent expansions — 
ultraviolet self-energies proportional to the length of the loops exponentiate and 
so are better handled by taking a logarithm.) Small Wilson loops are among 
the most ultraviolet quantities there are in lattice gauge theory. To probe lower 
momenta, we examine "Creutz ratios" of the loop values, 

, fW mn W m -i n -i\ . 

Xmn = - In — — . (41) 

\Wmn-iWm-inJ 

The most ultraviolet parts of the loops cancel in Creutz ratios and so these 
tend to be more infrared than the W's (ie, have lower <7*'s). In general quan- 
tities with larger loops and/or smaller g*'s are more likely to have significant 
nonperturbative contributions. 

2 This definition needs modification if there are nearly canceling contributions coming from 
very different scales. Luckily this situation does not arise often, and not at all in the applica- 
tions we are considering here. 
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lnWix oty(l/a) ay(n/a) 



a Ay 



5.4 
5.7 
6.0 
6.2 
6.4 
6.8 
9.0 
12 



0.7516 w 1 0.271 

0.5995 0.355 0.188 

0.5214 0.247 0.156 

0.4884 0.214 0.143 

0.4610 0.191 0.133 

0.4167 0.160 0.117 

0.2795 0.088 0.074 

0.1954 0.056 0.050 



4.1 x 10~ 3 

1.2 x 10~ 4 



0.627 
0.293 
0.169 
0.127 
0.097 
0.058 



Table 1: Monte Carlo data for logarithm of the plaquette, together with the 
coupling constant values and the ay scale parameter in lattice units. The inverse 
lattice spacing at (3 = 5.7 is a -1 « 1 GeV. All results are for the standard lattice 
action with n/ = 0. 

To make perturbative predictions we must determine ay for each value of (3 
used in the gluon action Eq. ( p8|) (ie, for each value of the lattice spacing). In 
principle we could do this by measuring the static-quark potential in the simu- 
lation. However it is simpler to redefine the coupling in terms of the expectation 
value of the plaquette operator. Using perturbation theory, one finds that 



where scale q* — 3.4/a is determined using the methods of the previous section, 
and n/ is the number of light quark flavors in the simulation (n/ = for the 
results we show here). This relation can be taken as a definition of ay. The 
coupling is determined nonperturbatively by inverting the series to determine 
ay(3.4/a) from simulated values for W\\. Having a value for ay(3.4/a), we 
obtain ay(q) at any other scale using two-loop perturbative evolution: 



where (3q = (11 - |n/)/4-7r, ft = (102 - 4^n/)/167r 2 , and Ay is a constant (the 
"scale parameter"). 

Simulation results for ay at several different fts are given in Table 0. Values 
for aAy are also listed. Since Ay is independent of ft these last numbers can 
be used to compute ratios of the lattice spacings at different fts. The inverse 
lattice spacing at (3 — 5.7 is about 1 GeV so that a » 0.2 fm. 

As a first test of perturbation theory, consider the Creutz ratio %22- It has 
the perturbative expansion 




V)}, (42) 



a v 1 (q) 



/3 ln(g7Ay) + /3i/ ) 9olnln( g 2 /A y ) + 0(ay) 
a v 1 (qo) + ((3o + Piav(q )) HlVlo) 



(43) 



X22 = 1-2 ay(l.l/a)- 0.4 o? v + 0{o? v ) 



(44) 
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av{^ /a) 


P = 5.4 


P = 5.7 


P = 6.0 


P = 6.2 


-InWn 


0.271 


0.188 


0.156 


0.143 


-I11W12 


0.269 


0.190 


0.156 


0.143 


- In W22 


0.271 


0.194 


0.156 


0.143 


-I11W33 


0.293 


0.205 


0.159 


0.145 



Table 2: Values of ay(ir/a) as determined by comparing three-loop perturbation 
theory with nonperturbative simulation results for the expectation values of 
small planar Wilson loops. 

where, as expected, its scale q* = 1.1 /a is significantly smaller than that for the 
plaquette. Monte Carlo simulation results are compared with first-order and 
second-order perturbative predictions for several /3's in Figure |l|. Perturbation 
theory gives excellent results throughout the range shown. 

For comparison results using a^gt^*) in place of ay(q*) are also shown; 
they are roughly comparable. Results are also shown for perturbation theory 
in terms of the bare coupling, «i at = gf at /4n. These are terrible. I show these 
last results only because it is a common practice in lattice perturbation theory 
to expand in terms of the bare coupling. This is a serious mistake — one that 
has in large measure crippled earlier attempts at using lattice perturbation 
theory. Using ai a t is wrong for two reasons. First, it ignores the fact that 
the Creutz ratio is significantly less ultraviolet than most Wilson loops and so 
should be expressed in terms of the running coupling evaluated at a smaller 
scale. Second, it ignores the renormalization in the relation between the lattice 
coupling and the continuum coupling (Eq. and, as we discuss in the 

next section, this renormalization is large. Expansions in ai at almost always 
underestimate perturbative effects. 

Our choice of the plaquette for determining ay is for convenience. In prin- 
ciple one ought to be able to use any short-distance quantity to measure ay in 
a simulation. If perturbation theory is working, the results obtained from dif- 
ferent quantities at a given P should agree. In Table | we list numerical values 
of ay(n/a) as determined from a variety of small Wilson loops at several values 
of p. In each case ay{q*) is adjusted so that the three-loop peturbative pre- 
diction agrees with simulation results, and then ay is evolved to scale q = it / a 
using the two-loop beta function. The agreement between different determi- 
nations of ay is good (even excellent for the smallest loops) for all /3's listed, 
including P = 5.4 where the lattice spacing is almost 1/2 fm. (Larger loops 
give poorer results because they are more strongly affected by nonperturbative 
effects.) 

The g-dependence of ay(q) can be determined directly from simulations. 
The data in Figure are obtained by fitting second-order expansions in ay(q*) 
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Figure 1 : Results for Creutz ratio X22 at different couplings (3 from Monte Carlo 
simulations (circles), and from perturbation theory (using ay{q*) (diamonds), 
a Ms(<7*) (boxes), and a\ at (crosses)). The first plot shows perturbation theory 
through one-loop order, and the second through two-loop order. Statistical 
errors in the Monte Carlo results are negligible. 
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to simulation results for the six smallest Creutz ratios and the six smallest 
Wilson loops. The value of ay{q*) obtained from each fit is plotted versus 
the q* for that quantity. The q*'s for the twelve quantities used here range 
from OA/a (for X44) to 3.4/a (for -lnWn). These simulation results may 
be compared with the two- loop perturbative prediction for ay(q) (solid line), 
arbitrarily normalized so that the curve passes through the point for — lnW22- 
The Monte Carlo data are quite consistent with perturbation theory, even at 
ft = 5.7. The smallest q*'s shown are on the order of several hundred MeV. 

These data, and much more (see later sections), show that lattice perturba- 
tion theory for QCD is very effective even at momenta as low as 1 GeV. Indeed 
it is about as effective as continuum perturbation theory. Perturbation theory 
seems to work so long as av{q) < 1. Possibly the most important consequence 
of this result is that it becomes likely that perturbative improvement of the 
QCD lagrangian will work even with lattice spacings as large as a = 1/2 fm. 
However, there remains one problem: perturbative expansions for the couplings 
in the corrected lagrangian tend to have very large coefficients. The origin of 
this problem and its cure are the subject of the next section. 

4 Tadpole Improvement 

The very large renormalization factor relating the bare lattice coupling a\ at to 
the continuum coupling ay is one of many examples where a large renormal- 
ization is required to relate a lattice quantity to its continuum analogue. These 
large renormalizations all have a common source: the compact nature of the 
link variable used in lattice theories. 

We design lattice operators by mapping them onto analogous operators in the 
continuum theory. For gauge fields, the mapping is based upon the expansion 

U^{x) = cvp(-igaA^) -> 1 - igaA^x). (45) 

This mapping seems plausible when a is small, but it is misleading in the quan- 
tum theory since further corrections do not vanish as powers of a. Instead one 
gets terms like ^{igaA^) 2 whose vacuum expectation value is suppressed only 
by g 2 — the explicit a 2 cancels a 1 /a 2 from the ultraviolet divergent (A 2 ) . These 
corrections turn out to be uncomfortably large. For example, 

<±TrE/ M ) LG w0.82 (46) 

at (3 = 5.7, while (| Tr (1 — igaA^)) = 1. These very divergent loop contribu- 
tions are referred to as "tadpole" contributions. 

The tadpoles spoil our intuition about the connection between lattice opera- 
tors and the continuum, leading to unexpected renormalizations. To regain this 
intuition we need to sharpen the relation between XJ^ and the continuum gauge 
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Figure 2: Values of ay(q) for a range of g's as determined from lattice QCD 
measurements at various (3's. The data points (circles) are measured values 
(with negligible statistical errors) obtained by fitting second-order perturbation 
theory to Monte Carlo simulation data for various short-distance quantities. 
The solid line shows the variation in oty (q) expected from two-loop perturbation 
theory. 
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field A^. Results like that in Eq. (46) suggest that the appropriate connection 
with the continuum is 



U,j.(x) — ► Mo (1 - igaA^x)) 



(47) 



where uo can be thought of as the mean value of the link operator (gauge invari- 
ance requires that it be an overall factor). This suggests that lattice operators 
will be more continuum-like if we replace every link operator by a "tadpole 
improved" operator [/„: 



but it is more convenient (and not very different) to define uq in terms of the 
expectation value of the plaquette: 



ing Uft about 15% smaller than U^. This difference is considerable, particularly 
in operators that involve products of several link operators. 

The tadpole-improved operators t/ M are much closer in their behavior to 
their continuum analogues than the unimproved operators; large tadpole renor- 
malizations are canceled by the uq. To see what effect tadpole improvement has 
on the design of lattice operators we now examine two applications. 

4.1 The gluon action 

Our new prescription for building continuum-like operators on the lattice sug- 
gests that 



with four powers of uq to correct for the four link operators, is a better action 
for lattice QCD. In particular, perturbation theory in ai a t = gi at /4w should 
be more like continuum perturbation theory (ie, no tadpoles). Of course, this 
action becomes the same as the standard action if we identify 



This relation is important because it explains why perturbative expansions in 
the bare coupling ai at fail so badly. Our tadpole analysis indicates that ai a t> 
not ai a t, is the correct expansion parameter for the theory. The difference is 
significant: for example, ai a t rs 1.8ai a t at (3 = 5.7. Using «i at almost guarantees 
that perturbation theory will give results that are much too small.tj 

3 An exception is the expectation value of the plaquette itself (without a logarithm). Be- 
cause this operator is identical to those in the lagrangian, tadpole corrections cancel in its 






(50) 



"lat = M aiat- 



(51) 
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If this analysis is correct, the tadpole-improved bare coupling aiat should be 
roughly equal to the continuum coupling avi^/a). We can check this in pertur- 
bation theory which gives a relation between ay and the usual bare coupling: 

a v {7r/a) = ^ {1 + A.7a v + 0(a 2 v )} . (52) 

We tadpole improve this expansion by dividing out the perturbative expansion 
for = 1/(| Tr[/Pjf q ) (Eq. (g§)) on the left-hand side: 

a ^/ fl ) = ,i T ti aq , {l + 0.513av + O(o^)} (53) 
= out {1 + 0.513a v + 0(a 2 v )} . (54) 

This equation implies that ay(-7r/a) and ai a t agree to within 10% at j3 = 5.7, 
while ai a t is almost a factor of two smaller. Thus perturbation theory confirms 
that almost all of the large renormalization factor relating ay to ai a t is due to 
tadpoles. 

Tadpole improvement is very important when improving lattice actions, since 
the correction terms tend to have lots of link operators. If you completed the Ex- 
ercises in the earlier sections, you found that the order a 2 errors in the standard 
lattice action for classical QCD are removed by replacing 

< aq - < aq - ^ K* 2 + Ulf) (classical), (55) 

where Uj^ 2 is the 1x2 planar Wilson loop. Since U^* 2 has two more link 
operators than i/P]f q , the tadpole improved correction in the quantum theory 
is given by 

^ aq - t^? q ~ ^| KZ 2 + Ulf) (tadpole-improved). (56) 

The extra factor of 1/uq is important, particularly at the large lattice spacings 
we want to use; without it the correction is too small and order a 2 errors are 
only partially removed. 

4.2 The quark action 

As a second example, consider the tadpole-improved version of Wilson's action 
for heavy quarks: 



X X,fl ^ ^ 



tp(x) + h.c. } . (57) 



expectation value. Thus the tadpole improved expansion and the expansion in terms of Qi a t 
are the same in lowest order, and naive, unrenormalized lattice perturbation theory gives 
pretty good results for this quantity. This is very much a special case. 
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Again, this action is identical to the usual one if we relate the modified param- 
eters, here the "hopping parameter" k, to the usual ones by rescaling with uq: 

k = K Uq. (58) 

Exercise: Setting U^/uq = 1, show that the quark mass in the free-quark theory is 
given by ma — 1 /2k — 4. 

The modified hopping parameter should be more continuum-like than the 
usual one; for example, the free-quark value that gives massless quarks, k c = 
1/8, should be roughly correct for interacting quarks as well. Thus a nonper- 
turbative formula for the critical value of the conventional hopping parameter 
is 

k c « l/8u 0l (59) 
which implies that the critical value of the bare quark mass is 

m c a = 1/2k c - 4 « 4(u - 1)- (60) 

Again we turn to perturbation theory to check the validity of this result. In 
perturbation theory, m c is given by 

m c a = -5.457a y (2.6/a) + 0(a v ) (61) 

where the scale of ay is determined as in Section |3.l| . Pulling out a 4(mq — 1) 
for the tadpole contribution, this becomes 

m c a = 4 ((iTrl^) 1 / 4 - l) - 1.268 a v (1.0/a) + 0{a 2 v ). (62) 

Thus perturbation theory confirms that the bulk of the renormalization of m c 
is due to tadpoles. Notice also that the scale in ay is smaller with tadpole 
improvement; this is because the very ultraviolet tadpole contributions have 
been removed. 

By removing the tadpole contributions, we make the perturbation theory 
more convergent, and therefore more accurate. The value of uq is obtained 
directly from the simulation (from the measured plaquette expectation value), 
and so we need never deal with tadpole effects in perturbation theory. The 
accuracy of this procedure is illustrated in Table |^ where we list values for 
m c a at different /3's. We give results from the ordinary perturbative expan- 
sion (Eq. (|6l|)), from the tadpole-improved expansion (Eq. ([32])), and from non- 
perturbative Monte Carlo simulations. The tadpole- improved results agree very 
well with the simulation results, and are about as accurate. 

Because there are link operators in the kinetic part of the quark action, tad- 
pole improvement of the links affects the relation between the lattice quark field 
and the continuum quark field. In the continuum limit, the tadpole-improved 
lagrangian for massless quarks becomes 

TK&t^ + O(o). (63) 
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a 


perturbative 


tadpole- improved 


simulation 


5.7 


-1.11 


-1.00 


-1.04(2) 


6.0 


-0.91 


-0.80 


-0.80(2) 


6.1 


-0.86 


-0.76 


-0.78(2) 


6.3 


-0.79 


-0.70 


-0.70(2) 



Table 3: Values for the critical mass m c a in Wilson's quark action as computed 
using ordinary perturbation theory, tadpole-improved perturbation theory, and 
Monte Carlo simulation. 

This indicates that the tadpole-improved lattice quark operator is 

ij) = \j2k c i\) = if}/ 2 (massless quarks), (64) 

where we use the fact that k c « 1/8. This lattice operator has roughly the 
same normalization as the continuum field; in particular, there are no large tad- 
pole contributions contributions to the renormalization constant relating them. 
This is important in designing new lattice operators involving quark fields. For 
example, if one wants to calculate matrix elements of the continuum current 
■07 At 7 5 'i/>, then one should simulate with the lattice operator 

■0 7 ^ 7 5 ^ = I^y^ 5 ^. (65) 

Exercise: Compute the tree-level correction needed to remove order a errors from 
the Wilson action. Tadpole-improve your result. 

5 Nonrelativistic QCD 

As an illustration of the ideas developed in these lectures we end with a review 
of recent simulation results concerning the upsilon family of bb mesons. This is 
an attractive system to study because so much is known about it. The quark 
potential model, for example, provides an accurate phcnomcnological model of 
the internal structure of the mesons. Also there is much experimental data for 
these mesons. The low-lying states are largely insensitive to light-quark vacuum 
polarization — for example, the T, T', \b- ■ -are all far below the threshold for 
decays into B mesons — and therefore can be accurately simulated with n/ = 0. 
Furthermore, these mesons are very small, the T being about five times smaller 
than a light hadron. This makes them an excellent testing ground for our ideas 
concerning large lattice spacings. 

The T spectrum indicates that the 6-quarks in the meson are nonrelativistic. 
This means that the most important momentum scales governing the meson's 
dynamics are smaller than the quark's mass Mb. We can take advantage of 
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this fact by choosing an inverse lattice spacing of order the quark mass, thereby 
excluding relativistic states from the theory. Then it is efficient to analyze 
the heavy-quark dynamics using a nonrelativistic lagrangian (NRQCD). The 
lagrangian used to generate the results shown below was 

£nrqcd = ^ (i-^)"^ (l-^) n (l-a«0^-Vty, ( 66 ) 
where n — 2, Hq is the nonrelativistic kinetic-energy operator, 

A( 2 ) 

Ho = (67) 

M 6 ° is the bare quark mass, and SH is the leading relativistic and finite-lattice- 
spacing correction, 



8(M°) 3 \ 2n J 24M° 

• B + (A • E - E • A) 



2M° 8(M°) 2 

-8(]|F <r - (AxE - ExA) - (68) 

Here A and A^ 2 ) are the simple gauge-covariant lattice derivative and laplacian, 
while A^ 4 ) is a lattice version of the continuum operator JZ-Df. The chromo- 
electric and chromomagnetic fields, E and B, are defined in terms of different 
components of 

U$«-U$«\ (69) 

The entire action was tadpole improved by dividing every link operator 
by uq. Potential models indicate that corrections beyond SH contribute only 
of order 5-10 MeV to T energies. The simulation results shown below used 
gluon fields generated at j3 = 6, which corresponds to a lattice spacing of about 
1/12 fm or about half the radius of an T. 

The details of this lagrangian are unimportant to us here. What matters is 
that SH consists of correction terms just like the ones we have been analyzing 
for other theories, the only difference here being that we are correcting both for 
finite-a and for the absence of relativity (ie, order v 2 /c 2 errors). The coefficients 
of the correction terms were determined with tree-level perturbation theory 
and tadpole improvement, using precisely the techniques outlined in the earlier 
sections. So the extent to which SH improves the simulation results is a measure 
of the efficacy of all the techniques discussed in these lectures. 

The T spectrum is very well described by the simulation. Simulation results 
for the low- lying 3 Si and 1 P\ energies are shown in Figure |. These compare well 
with experimental results (the horizontal lines), as they should since systematic 
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Figure 3: NRQCD simulation results for the spectrum of the T( 3 Si) and /ib( 1 Pi) 
and their radial excitations. Experimental values (dashed lines) are indicated 
for the 5-states, and for the spin-average of the P-states. The energy zero for 
the simulation results is adjusted to give the correct mass to the T. 



errors are estimated to be less than 20-40 MeV. It is important to realize that 
these are calculations from first principles. The only inputs are the lagrangians 
describing gluon and quark dynamics, and the only parameters are the bare 
coupling constant and the bare quark mass. In particular, these results are not 
based on a phenomenological quark-potential model. These are among the most 
accurate lattice results to date. 

The corrections in 5H have only a small effect on the overall spectrum, 
but the spin structure is strongly affected. The lagrangian without 5H is spin 
independent, and gives no spin splittings at all. Simulation results for the spin 
structure of the lowest lying P state are shown in Figure [|. Again these compare 
very well with the data, giving strong evidence that corrected lagrangians work. 
Systematic errors here are estimated to be of order 5 MeV. Note that the spin 
terms in 5H all involve either chromoelectric or chromomagnetic field operators. 
These operators are built from products of four link operators and so tadpole 
improvement increases their magnitude by almost a factor of two at /3 = 6. 
Without tadpole improvement of the spin splittings from this simulation would 
have been much too small. 

The correction terms also affect the simulation mass of the T. Our non- 
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Figure 4: Simulation results for the spin structure of the lowest lying P-state 
in the T family. The dashed lines are the experimental values for the triplet 
states, and the experimental spin average of all states for the singlet (hb). 



relativistic action gives only part of the total energy or mass of the meson. 
The full mass is obtained by adding the masses of the quarks to the meson's 
nonrelativistic energy -Enr : 

M t - 2 (Z m Ml - E ) + £ N r (70) 

where Z m and Eq are ultraviolet renormalizations that are computed using 
perturbation theory, M® is a tunable parameter of the theory, and -Enr is com- 
puted in the simulation. The bare quark mass for the simulation discussed 
here was tuned so that this formula gave Mr = 9.5(1) GeV (the correct value 
is 9.46 GeV). The T mass was also determined a second way in the simulation by 
computing the nonrelativistic energy of an T as a function of its momentum p, 
and fitting its low-momentum behavior to the form 

E r (p) = E NR + -^— +-... (71) 

This determined the kinetic mass M^ n of the meson. In a purely nonrelativis- 
tic theory the kinetic mass equals the sum of the quark masses. Only when 
relativistic corrections are included is this mass shifted to include the binding 
energy, giving Mr- The simulation without 5H gave Mki n = 8.2(1) GeV, which 
is quite different from the upsilon mass determined using Eq. (|7C|). With SH, 
the simulation gave Mki n = 9.5(1) GeV, which is in excellent agreement. All of 
the spin- independent pieces of SH contribute to the shift in Mk; n ; once again 
we have striking evidence that corrected actions work. 

This simulation has only two parameters: the bare coupling constant and 
the bare quark mass. These were tuned to fit experimental data. From the 
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bare parameters we can compute the renormalized coupling and mass. This 
simulation implies that the renormalized or "pole mass" of the 6-quark is 



M b = 4.94(15) GeV. 



(72) 



The renormalized coupling that is obtained corresponds to 



a&UMz) = 0.112(5), 



(73) 



which agrees with results from high-energy phenomenology and is about as ac- 
curate. This last result is striking: it shows that the QCD of hadronic structure 
and the QCD of high-energy quark and gluon jets are really the same theory. 

6 Conclusion 

The T simulations described in the last section show that lattice QCD can 
produce accurate results even when the lattice spacing is as large as half the 
radius of a meson. For light hadrons this would correspond to a w 0.4 fm or 
j3 w 5.4 (for simulations with n/ = 0). Our analysis shows that perturbation 
theory is still reliable at such distances, and so perturbative improvement of 
the lagrangians used in simulations probably still works. Given that finite- 
volume errors become manageable for lattice sizes of order 2-3 fm, it seems 
likely that reliable simulations of full QCD are possible on lattices as small 
as 6 4 . Simulations on such small lattices are literally a thousand times faster 
than simulations on the 20 4 lattices commonly used today. If coarse lattices 
really do work, the shift to small lattices and improved lagrangians will have a 
revolutionary effect on numerical QCD. 
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